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The entropy properties are analyzed by Ma's coincidence method in 
7r + p and K + p collisions of the NA22 experiment at 250 GeV/c incident 
momentum. By using the Renyi entropies, we test the scaling law and 
additivity properties in rapidity space. The behavior of the Renyi entropies 
as a function of the average number of particles is investigated. The results 
are compared with those from the Pythia Monte Carlo event generator. 

PACS numbers: 13.85.-t, 13.85.Hd, 05.10.-a 

1. Introduction 

It was suggested recently that the event coincidence probability method 
of measuring entropy, originally proposed by Ma is well suited for the 
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analysis of local properties in multiparticle systems produced in high energy 
collisions [21 El EI 03 El El El ■ In heavy ion collisions it can be used in the 
search for formation of a quark gluon plasma (QGP) [2], but the method is 
equally applicable to hadron-hadron collisions It also proved effective 
for the study of systems of particles created in multiplicative branching 
processes as present in high-energy QCD jet fragmentation. Finally, 
entropy measurements offer an additional tool to analyze event-to-event 
fluctuations and particle correlations [U E] ■ 

The existence of a QGP at high energy density is a prediction of QCD. 
Many theorists and experimentalists are searching for possible signals of the 
formation of such a new state of matter. The primary goal of the Relativistic 
Heavy Ion Collider (RHIC) at Brookhaven National Laboratory is to create 
and to study this deconfined state. Systematic measurements of the local 
entropy at RHIC may provide direct information about the internal degrees 
of freedom of the QGP state and its evolution. It is thus very important 
to first see how the proposed entropy measures behave in hadron-hadron 
collisions. 

In this note, we present the results of a study of entropy in 7r + p and K + p 
collisions at 250 GeV/c incident momentum [yfs = 22 GeV) collected by the 
NA22 experiment. We investigate the dependence on discretization of the 
system, the effects of particle correlations by testing scaling and additivity 
properties of the entropy measures, as well as the multiplicity dependence 
of the Renyi entropies in rapidity space. We also compare the results of the 
NA22 data to those of the Monte Carlo event generator Pythia 5.7 used 
for inelastic, non-single-diffractive pion-proton collisions with Bose-Einstein 
Correlations (BEC) included, as well as to a simple random production 
model. We end with a short conclusion and outlook. 



The original definitions of the standard Shannon entropy and the Renyi 
entropies are: 



Here, pj = nj/N denotes the probability to obtain a specific configuration of 
the system, where rij is the number of events in such a configuration and N 
is the total number of events. The sum runs over all possible configurations. 

Following Ma's method [Q and as explained in j^j, the first step in the 
entropy measurement is to determine coincidence probabilities. For every 
event, a certain phase space region is divided into M bins of equal size. An 
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event is then characterized by the number of particles, m,;, in each bin, i.e., 
by a set of integer numbers s = {mi}, where i = 1, ...M. 

After counting how many times, n s , the set s appears in the whole event 
sample, one can determine the sums 

N k = J2n s (n s -l)...(n s -k + l). (2) 

s 

This gives the total numbers of observed coincidences of k configurations. 
The coincidence probability of k configurations is then given by 

° k = N(N -l)...(N -k + 1) ' ® 

where N, as above, is the total number of events in the sample. Only states 
with n s > k contribute to C k . 

From the coincidence probabilities, one calculates the Renyi entropies 
defined as JO] 

The Shannon entropy S is formally equal to the limit of H k as k — > 1. 
So it can only be obtained by an extrapolation method. One possibility we 
will investigate in some detail is to take [3] 

H k = a-j^j + °o + a x {k - 1) + a 2 {k - if + (5) 

These expansions may be used as a system of linear equations to deter- 
mine a and cij, with as many terms as the number of H k s available. The 
Shannon entropy is then calculated as S = a + ao- It is instructive to see if 
the value of S depends on the number of Renyi entropies used. In general, 
it turns out that already the value obtained from two terms, S[H2, H3], and 
that from three terms, S[H%, H3, H4], nearly coincide. 

One of the most attractive features of Ma's coincidence method com- 
pared to the standard method of using Eq. (1) directly is that, as seen from 
Eq. (3), the statistical error decreases very fast with increasing number of 
configurations. Moreover, it has been demonstrated jH] that the extrapola- 
tion method yields more stable results than the standard method. 

For a system close to thermal equilibrium, and if the phase space sub- 
division is sufficiently fine-grained, the scaling relation 



H k (lM) = H k (M) + InZ =>- S(IM) = S{M) + In/ 
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holds. Here, M and IM are the numbers of bins in two different discretiza- 
tions. 

Another feature is additivity: the entropies measured in a region R 
which is the union of two non-overlapping and independent regions R\ and 
R 2 satisfy 

Hk{R) = Hk(Ri) + Hk(R 2 ) 

=> S(R) = S(R 1 ) + S(R 2 ). (7) 

We will check these two important properties below. 

3. Data sample 

The NA22 experiment made use of the European Hybrid Spectrometer 
(EHS) in combination with the Rapid Cycling Bubble Chamber (RCBC). 
Details on the setup and on the reconstruction of the data can be found 
in Q31 1121- 
Charged-particle momenta are measured over the full solid angle with 
an average resolution varying from 1-2% for tracks reconstructed in RCBC 
and 1-2.5% for tracks reconstructed in the first lever arm, to 1.5% for tracks 
reconstructed in the full spectrometer. Ionization information is used to 
identify and exclude protons up to 1.2 GeV/c and electrons (positrons) up 
to 200 MeV/c. All unidentified tracks are given the pion mass. 

In our analysis, events are accepted when the measured and recon- 
structed charged-particle multiplicity are the same, no electron is detected 
among the secondary tracks and the number of badly reconstructed tracks 
is 0. The selection of the data is described in detail in ^Hj. After all nec- 
essary rejections, a total of 44,524 inelastic, non-single-diffractive events is 
obtained. Acceptance losses are corrected by a multiplicity-dependent event 
weighting procedure. 

4. Results and Discussions 

In Fig. 1, we show the Shannon entropy and the Renyi entropies H 2 to 
Hi as obtained from an M = 9-fold division of central cms rapidity windows 
of size Ay ranging from 1 to 6 units. The Shannon entropy is calculated 
by extrapolation according to Eq. (5). The solid symbols show the NA22 
data, the open symbols the Pythia predictions. All Renyi and Shannon 
entropies increase as the rapidity window widens. The values of Hf. decrease 
with increasing k. The figure further shows that the two-term extrapolated 
Shannon entropy S[H 2 , H3] agrees well with the three-term extrapolation 
S[H 2 , H$, #4]. For the Hk, Pythia tends to be flatter than the data. This 
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Fig. 1. Shannon and Renyi entropies obtained from an M = 9-fold division of 
central rapidity windows of size Ay ranging from 1 to 6 units. 



trend increases with increasing order k, so that the extrapolation procedure 
leads to a Shannon entropy S which rises more steeply than the data. 

It may be worthwhile to point out that, in general, the extrapolation 
procedure for the Shannon entropy is not unique and thus will introduce an 
additional uncertainty |HJ HJ. For example, using a polynomial expansion 
(Eq. (20) of ref. [3]) for Ay = 6 results in a value which is about 10% 
smaller for S[H2, H3] and 7% smaller for S[H2, H^,H/^\. The Renyi entropies 
are of great interest by themselves and provide valuable information about 
the system without the need for any extrapolation. In the following, we 
therefore concentrate on the Renyi entropies, in particular on H2. 

If scaling according to Eq. (6) holds, one should observe a linear relation 
if Hk(M) is plotted as a function of — InJy, where 5y = Ay/M is the bin 
size. 

We first test scaling of H2 in central rapidity windows in Fig. 2. The solid 
circles, squares and triangles represent the NA22 data for \y\ < 0.5, 1.0, 2.0, 
respectively. They all are flattening with increasing — InJy, i.e., the scaling 
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Fig. 2. The dependence on the dividing bin sizes of the second Renyi entropy 
calculated in central rapidity windows. 



law does not hold for the NA22 data. This implies that there are strong 
particle correlations in the process under investigation. The open symbols 
in Fig. 2 represent the Pythia results with a Bose-Einstein correlation 
strength A = 1.0 9 . They do not fully agree with the data. We have 
verified that this discrepancy is not sensitive to the actual value of A used. 

The cross symbols for \y\ < 1.0 in Fig. 2 show the results from a random 
model generated as follows: take the multiplicity distribution in the window 
\y\ < 1.0 from the data, but the multiplicity fluctuations in each 5y bin 
to be Poissonian with a mean as in the data. In this model, there are 
no correlations between particles, so it should give a nearly straight-line 
relationship, as indeed shown by the figure. 

Due to the asymmetry of the proton and meson fragmentation regions in 
the NA22 experiment, we plot forward and backward hemispheres separately 
in different non-central windows of fixed size Ay = 1 in Fig. 3 (a) and (b), 
respectively. We see that, again, the H2 values are flattening with increasing 
— ln<5y, so that no linear relationship is found. Pythia has the same trend, 
but overestimates the data in the peripheral regions. We also note that H2 
becomes smaller for more peripheral windows, due to decreasing particle 
density, i.e., decreasing coincidence probability. 

In order to test additivity, Fig. 4 (a) shows H2 obtained in two ways: one 
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Fig. 3. The dependence on the dividing bin sizes of the second Renyi entropy 
calculated in (a) forward and (b) backward hemispheres in different non-central 
rapidity windows. 
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Fig. 4. The second Renyi entropy as a function of — In 5y measured in two (a) 
adjacent and (b) separate intervals. 



is a direct measurement in the rapidity window \y\ < 1 (solid circles), the 
other the sum of the values obtained from the adjacent ranges — 1 < y < 
and < y < 1 (solid squares). A clear difference is observed between 
these two results. The H2 obtained from summation is larger than the one 
calculated directly, and the difference increases with increasing — InSy. So, 
additivity of H2 is not observed here. The reason is that strong correlations 
exist between the particles belonging to these two adjacent y ranges. Indeed, 
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Table 1. Fit values of Renyi entropies in different rapidity windows. 
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1.79/8 








Hz 
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0.50±0.03 


4.22/8 



no violation of additivity is observed for the random model (crosses and 
diamonds in Fig. 4 (a)), thus confirming the role of correlations in the data. 

Fig. 4 (b) shows the results from two regions separated by two rapidity 
units. In this case, the difference between the direct measurement in the 
window 1 < \y\ < 2 (solid circles) and the sum of the two non-adjacent 
windows —2 < y < — 1 and 1 < y < 2 (solid squares) is very small. PYTHIA 
(open symbols) shows the same trend but again overestimates the values in 
the non-central rapidity regions. 

The additivity observed here is consistent with the experimental ob- 
servation that correlations between particles belonging to widely separated 
rapidity intervals are small [Tl|- However, at much higher energies, strong 
long-range correlations are observed ^] which should lead to a breaking of 
the additivity property. This was confirmed in Pythia model calculations 
for proton-proton collision at yfs = 1800 GeV 

As advocated in the dependence of Renyi entropies on particle multi- 
plicity n carries important information on the produced system. An entropy 
proportional to n is indicative of an equilibrated system with no strong long- 
range correlations. On the other hand, for non-equilibrated systems or a 
superposition of sub-systems with different properties, proportionality to 
In n is expected. 

In Fig. 5, we show the charged-particle multiplicity dependence of H 2 
(solid circles) and H 3 (solid squares) for three different rapidity windows, 
where the number of divisions M is chosen such that 5y is the same for each 
rapidity window. Table I shows the fit results for the logarithmic function 

V = ci +c 2 lnn ch . (8) 

In the central region, a logarithmic rather than a linear relation is observed. 
This again confirms that there is no thermal equilibrium in the system. 
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Fig. 5. The dependence of the second and third Renyi entropies on multiplicity in 
linear (left) and logarithmic (right) scale, with (a)(b) M = 6-fold division in the 
rapidity window \y\ < 3, (c)(d) M = 4 in \y\ < 2 and (e)(f) M = 2 in \y\ < 1. 

Pythia (open symbols) agrees quite well with the data. 



5. Conclusions and outlook 

We have analyzed the entropy properties of multiparticle production 
in 7r + p and K + p collisions at 250 GeV/c. By using Renyi entropies, we 
find that both the scaling and the additivity properties are not generally 
valid and that the multiplicity dependence shows a logarithmic rather than a 
linear relation. All of these confirm the expectation that thermal equilibrium 



10 



ENTROPY PRINTED ON FEBRUARY 7, 2008 



is not reached in hadron-hadron collisions at yfs = 22 GeV. 

The Pythia Monte Carlo model agrees, in general, quite well with the 
data. However, significant deviations exist. In particular, the model overes- 
timates the values in the peripheral rapidity regions, presumably due to too 
weak correlations. This shows that Renyi entropies provide a new sensitive 
measure of multiparticle correlations. 

It would be interesting to investigate the entropy properties of high- 
temperature, high-density systems, which may create the long expected 
QGP. RHIC has already collected data on gold-gold collisions at 130 and 
200 GeV and pp collisions at 200 GeV. The results presented here should 
provide a valuable guide to the interpretation of these measurements. 
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